Onsager relations in a two-dimensional electron gas with spin-orbit coupling 
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Theory predicts for the two-dimensional electrons gas with only Rashba spin-orbit interaction a 
vanishing spin Hall conductivity and at the same time a finite inverse spin Hall effect. We show 
how these seemingly contradictory results are compatible with the Onsager relations: the latter do 
hold for spin and particle (charge) currents in the two-dimensional electron gas, although (i) their 
form depends on the experimental setup and (ii) a vanishing bulk spin Hall conductivity does not 
necessarily imply a vanishing spin Hall effect. We also discuss the situation in which extrinsic spin 
orbit from impurities is present and the bulk spin Hall conductivity can be different from zero. 

PACS numbers: 72.25. Ba, 72.25 Dc, 72.25.Rb 
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R has been repeatedly questioned in the literature 
whether the Onsager relations [l| between direct and in- 
verse spin Hall effect are satisfied. In particular it has 
been argued that with the conventional definition of a 
spin current — defined as the product of spin and ve- 
locity operators — one cannot establish an Onsager rela- 
tion 0. Most recently doubts about their validity have 
been formulated Q after the prediction of a finite in- 
verse spin Hall effect in the two-dimensional electron gas 
with Rashba spin-orbit (SO) coupling, Ijl a system where 
the spin Hall conductivity vanishes |5l4lO|. In this paper 
we will cast the SO interaction in terms of non-Abelian 
SU{2) gauge fields [lll - flBj and show that: i) Onsager 
relations do hold in the presence of SO coupling, pro- 
vided the appropriate form of the spin current is used 
- crucially, this will depend on the particular measuring 
scheme employed; ii) a vanishing bulk spin Hall conduc- 
tivity does not imply a vanishing spin Hall or inverse 
spin Hall effect. We will discuss in some detail the ex- 
perimental relevance of our results, which will be shown 
to be valid in the presence of extrinsic SO coupling from 
impurities, too. 

To begin our discussion, let us imagine a two- 
dimensional electron gas (2DEG) with SO coupling. The 
Hamiltonian is 



add a time dependent perturbation Vi(<) of the form 



(1) 



imp is a random potential due to impurities. 



where 

taken to be s-wave scatterers. For definiteness' sake 
we choose for Hso the Rashba SO interaction, Hso = 
—apxT^ + apyT^ , though any other linear-in-momentum 
SO term could be handled (see below); and are 
Pauli matrices, a is the SO coupling constant. We now 
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eA,{t) + —TiA^it) 
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The vector potential Ai{t) is related to the electric field 
via Ei = —dfAi and is coupled to the particle current 
ji, whereas Al{i) is a fictitious spin dependent vector 
potential which creates a spin electric field = —dtAf 
and which is coupled to the conventional spin current jf ; 
77 is a formal SU{2) coupling constant. Physical mech- 
anisms actually generating this type of spin-dependent 
vector potential are discussed in Ref. [^. Lower [upper] 
indices indicate real space [SU(2), i.e. spin space] com- 
ponents. The Onsager relations connect the spin current 
generated by an electric field to the particle current gen- 
erated by a spin-electric field. For the spin Hall effect we 
conclude from Eq. ^ 

(4) 

where a^^ is the spin Hall conductivity and lu is the 
frequency. Instead of introducing the electric field E via 
the vector potential A(<), one could, cquivalently, choose 
a scalar potential 0(x, t) ~ — E(t) -x. One could then ask: 
will the Onsager relations ^ still hold once the spin- 
electric field is introduced via a spin dependent scalar 
potential? With the conventional definition for the spin 
current introduced in Eq. ([3]) the answer is "yes" only 
for vanishing SO coupling. This means that for a = 
the time dependent perturbation 
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with the spin dependent scalar potential — —E^ ■ x 
generates the same currents as V\. Formally, the two 
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cases Hi = H + Vi, H2 = H + V2 are connected by a 
gauge transformation 







Hi = UH2W 



ihUdtU^ 



(6) 



with U — exp[— (ir^/2/i)7yx(x, t)] and (9fx(x, t) = 
\I'^(x, <). On the other hand when a ^ Hi and H2 
are not connected by any gauge transformation. This 
can best be seen by writing the SO couphng in terms of 
a spin-dependent vector potential, 



(7) 



where the subscript R stands for "Rashba" and ri{AR)rl. = 



'2m.a, ri{AR)l 



2ma. Notice that within this ap- 



proach a different SO interaction - e.g. Drcsselhaus, a 
spatially modulated Rashba and so on - could be treated 
just the same and would simply amount to a different 
choice of SU{2) gauge fields. Now the external fields 
Vi and V2 are not equivalent any more, since under the 
gauge transformation ^ 



a^Q H + V2 ^ H' + Vi ^ H + Vi, 



(8) 



i.e. U sends V2 — > Vi and at the same time rotates the 
background Rashba field tjAr sending H H' . Explic- 
itly, to first order in x the spin dependent vector potential 
changes as 

vA'U^,t) = ^A'i + vfi-~\Me^'^A\ ~ <5°^77a,,x(x,t), 

(9) 

where e°'^^ and 5°'^' are the fully antisymmetric Ricci ten- 
sor and the Kronecker delta. The Rashba SO term is 
modified due to the second term on the right-hand side 
of Eq. (|9]). Physically this is unacceptable: the back- 
ground Rashba field has to remain fixed, else we would 
be describing a different system. Such a change can how- 
ever be absorbed into a redefinition of the spin current: 
fixing the background vector potential Ar, requires us to 
modify the definition of the current coupled to the exter- 
nal perturbation. To appreciate this point let us take 



(10) 



and gauge transform H2 using U previously defined. To 
linear order in the spin-electric field the result is 



H + V2{t)^H + Vi{t), 



with 



Vl(t) 



it 



triEl 



(11) 



(12) 



where jT is the conserved spin current operator suggested 
in Ref. |2|. Reintroducing the U{1) electric field we can 
write the equivalent of Eq. ([3]) 



Vlit) = J2[h^Mt)+J-r]At{t)] 



(13) 



and immediately obtain the Onsager relations 



-sH 



{uj)eEy{uj). 



(14) 

Eqs. Q and are the first main result of this work. 
They show that Onsager relations do hold in the presence 
of spin-orbit coupling, but the quantity reciprocal to the 
particle current changes depending on the experimental 
setup - i.e. on the way the external spin-electric field is 
generated. This means that the transport coefficient, the 
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spin Hall conductivity, changes too 

For lincar-in-momentum SO interaction the specific 
form of the spin Hall conductivity can be computed for 
any kind of spin electric field relying on the microscopic 
formalism developed in Ref. [l^ . which we will now fol- 
low. The goal is to verify explicitely the Onsager rela- 
tions Eqs. (HI and (HH). Let us then focus on the diffu- 
sive regime, in which the equations acquire a remarkable 
physical transparency. Generally, the particle and spin 
currents are the sum of a diffusion, drift, and a Hall cur- 
rent, the latter being responsible for the Hall and spin 
Hall effects. For a system without inversion symmetry, 
as it is the case for the Rashba model, extra terms ap- 
pear since a homogeneous non-equilibrium spin density 
can generate a spin current. In the SU{2) formulation 
such extra terms are automatically built in, and the par- 
ticle and spin current densities read [l^ 



j = -i?Vp + aE--L Vr xB^ (15) 

m ^ — ^ 

(16) 
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when the conventional definition of the spin current is 
used. Here, D = VpT/2 is the diffusion coefficient, Nq 
the density of states at the Fermi level, r the elastic scat- 
tering time and a = —2eNoD, i.e. the electrical conduc- 
tivity up to a charge — e. The above equations have been 
derived under the assumptions of weak disorder epT ^ h 
and weak SO coupling app f^/T- and pp being the 
Fermi energy and momentum, respectively In the fol- 
lowing, for simplicity, we will use units such that h = \. 
The SU (2) nature is manifest in the covariant derivative 
[Vis]° = dxiS'^ — e'^'"^r]A'^s'^ and in the spin dependent 
electric and magnetic fields 



- -dtA'} - a,.*" - e^^^r^^M^, 



^abc 



VA'^AI) 



(17) 
(18) 



For the Rashba model there is only one nonvanishing 
field, namely rjB^ = —{2'ma)^. Adding the external per- 
turbations Vi or V2 introduces further fields. We first 
consider Vi, Eq. ([2]), and obtain the additional fields as 
f| = iojA^, B% ~ —{2ma)A%, having moved to Fourier 
space {dt — icj, V — > iq) for later convenience. In 
linear response to the perturbation Vi, the transverse 
particle current generated by the spin-electric field £^ is 



3 



(about this point we disagree with Ref. 0, see also the 
appendix of this paper) 



7 =^B^7^ 



(19) 



where the dimensionless number 7 ~ —ma^T = 7i,it char- 
acterizes the couphng strength between spin and particle 
currents. A non zero spin-charge coupling signals the 
occurrence of the spin Hall effect [l3| independently of 
the spin Hall conductivity being different from zero or 
not, the latter fact depending of the experimental setup 
and other possible interactions in the Hamiltonian. The 
expression for the spin current of Eq. (jl6p reads 



(20) 



and in order to find its value we need the spin densi- 
ties. These can be obtained by solving the associated 
diffusion equations, which are nothing but the continu- 
ity equations for the currents (fTSj) an d_([l6|) . provided the 
SU{2) covariant derivatives are used[14| 



(21) 



with [dtsY = 9(5° -I- e°'^'^r]'^'^ s'^ . In particular the equa- 
tions for the in-plane spin densities in Fourier space are 
(^' = for the present case of Vi) 



iq • + 2majl = 
iq • + 2majy = 0. 



(22) 
(23) 

Inserting the Fourier transform of Eqs. p ^ -(fT6 |) into 
(P ^ - (|23p one obtains in the spatially homogeneous sit- 
uation 



4e 



(24) 



'DP 



where we have introduced the Dyakonov-Percl spin re- 
laxation time T^p = {2ma)^D. Notice that Eq. is 
non analytic in w and r^p . In the absence of Rashba SO 
coupling, i.e. in the limit r^p — >■ 0, the spin current is 
given by the spin electric field according to Ohm's law. 
When SO coupling is present, the spin current vanishes in 
the DC limit, i.e. w — !• 0. In the appendix this is shown 
explicitly by evaluating the Kubo formula diagrammat- 
ically. Relation yields the particle-current response 
to the spin-electric field and, to leading order in , 
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(25) 
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As required by the Onsager relations (|4]) this agrees with 
the spin Hall conductivity determined by the response 
of the conventionally defined spin current to the electric 
field. The latter result can be obtained by combining the 
expression for the spin current (jl6|) with the continuity 
equation 



We can now follow the same route while consider- 
ing the external perturbation V2, Eq. ([5|), with an x- 
dependent spin-scalar potential ^^{x,t). The latter in- 
troduces the following fields = £jj = —2ma'^^, f| = 
^iqx"^^- Our system is now homogeneous only along y, 
and the diffusion equations read 



js- = {-Dql - T^\,)s- + 4mal?zfe[s- - {Nol2)i^^'], 



'DP 



2r^'p)[s^ - {No/2)r^^' 



(26) 

AmaDiqxS^ , 



where we have ignored all terms that are quadratic in 
the external field notice that in the absence of no 
spin polarization exists, thus the spin density is itself at 
least 0{'^^). Solving Eqs. ([26|) for a homogeneous but 
frequency dependent spin-electric field we find 



cr?/ 



'DP c-z 
-1 "-"XJ 
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2tjjp —iuj 



(27) 



'DP 



and with Eq. 
ductivity is 



T9|) we conclude that the spin Hall con- 
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'DP 



e —iuj + 2Tpip —iuj - 



(28) 



'DP 



According to the Onsager relations (|T4)) . the reciprocal 
quantity to the inverse spin Hall current jy is the con- 
served spin current j| generated by an homogeneous and 
frequency dependent electric field along y. 



Jx 



= hm — s 

<?i->0 q^ 



(29) 



The diffusion equations to solve are now [we drop terms 



tuts = — T 



DP 



+ 2ma + "fcrEy 



(30) 



-iujs^- = -2T^pS^ - iqx [2{2ma)Ds^ — "faEy] . 

Their solution yields = a^^ {ijj)eEy with a'^^ [u) given 
by ([25]), thus verifying the validity of Eq. ([Til). 

It is now worthwhile investigating the robustness of 
the above results to the presence of extrinsic SO inter- 
action arising from impurities, since the latter are usu- 
ally present in real samples and in this case the static 
spin Hall conductivity {uj — > 0) is different from 
zero [3 [1^. To this end we add to the Hamiltonian the 
extrinsic term 



)2 

^(T X VT^mp 



W • p, 



(31) 



where Ao is the effective Compton wavelength describing 
the SO coupling in the system. The extrinsic SO inter- 
action pi[) modifies the theory only in two main aspects. 
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First, the presence of the extrinsic SO scattering intro- 
duces the Elliott- Yafet spin relaxation time, r,, so that 
Eq. ^ is modified to 



(32) 



with Ts — t(AoPf/2)^^. The second ingredient is that 
the parameter 7 entering Eq. ([TO)) acquires a contribution 
from the skew-scattering and side-jump mechanisms 7 = 
7int + 7ss + Jsj where 7int = —ma^T as before, while 
Isj = {\o/2f{m/T) and 7^5 = -(AoPF/4)^(27riVo?ii), ■!;o 
being the impurity scattering amplitude, see Ref. [20[ for 
technical details. One can now proceed as before and 
check that in linear response to Vi and V2 the relations 
(HI and (fT4)) still hold, with the spin Hall conductivities 
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-iuj + T- 



e -iiu + Tjjp + 



-1 



70" 



-iuj 



(33) 
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^' DP 



-IUJ - 



-■ (34) 
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We wish to stress two important points. First, in ob- 
taining the above we could still exploit Eq. ([29| . since 
the spin current j° introduced in Eq. (IT^ is by defini- 
tion conserved with respect to the full background field 



Hso + -f^cxtr[2l|. Second, and experimentally important. 



in the absence of intrinsic SO coupling, one has to take 
the a — > limit first, so that {uj) = a-*^(w), i.e. 
the two experimental setups corresponsing to Vi and V2 
become equivalent, since the out-of-plane spin density 
becomes a conserved quantity [26[. This is not the case 
in the presence of both intrinsic and extrinsic SO mech- 
anisms, since Vi is capable of sustaining a steady state 
bulk spin Hall current, whereas V2 is not. 

The relevance of our results with respect to available 
experiments is worth a more detailed discussion. Theory 
tells that in the pure Rashba case the bulk spin Hall con- 
ductivity vanishes, it is neither possible to drive a spin 
current by a uniform and weakly time dependent elec- 
tric field, nor to drive a charge current by (i) a uniform 
but weakly time dependent spin-vector potential (ii) a 
weakly space dependent but static spin-scalar potential. 
On the other hand when both intrinsic and extrinsic SO 
interaction are present, the bulk spin Hall conductivity 
can be different from zero. A cubic Dresselhaus term has 
a similar effect 0, [l^]. However even a vanishing bulk 
spin Hall conductivity does not imply the absence of the 
spin Hall effect and its inverse. The spin Hall effect and 
an induced edge spin-polarization are present close to a 
interface where non-spin-polarized carriers are injected 
into the Rashba 2DEG. This has been predicted first in 
Ref. and verified numerically in Ref. 23. This is also 
manifest in the expression for the spin current, Eq. (|16p . 
since when spin polarization is negligibly small the cur- 
rent becomes 



4™ 



(35) 



For the inverse spin Hall effect the situation is analo- 
gous. In an experiment such as the one of Ref. |2J no 
spin-electric field is applied to the samples. Instead a 
circularly polarized laser beam is used to create electron- 
hole pairs at a p-n junction between a 2DEG and a two- 
dimensional hole gas. With the junction suitably biased, 
spin polarized electrons are injected in the 2DEG, so the 
spin current j° at the interface is directly determined by 
the experimental setup and thus creates a Hall signal. 



jHall 



m ^ — ' 



(36) 



To conclude, we have shown the existence of Onsager 
relations connecting electric to spin-electric stimuli in a 
two-dimensional electron gas with spin-orbit coupling. In 
order to be explicit we focused on the Rashba model, but 
the non-Abelian formulation employed can be used for 
any linear-in-momentum SO interaction, possibly slowly 
varying in time and space, too. Quite important from 
the experimental point of view, the Onsager relations ob- 
tained are robust to the inclusion of extrinsic SO coupling 
from impurities and their specific form depends crucially 
on the measuring scheme employed. 
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APPENDIX - In the main text we have shown that a 
static spin-electric field introduced via a perturbation Vi 
does not create a spin current, a result which docs not 
agree with Eq. (15) of Ref. 0. To further support this 
statement we show here how to obtain this result with 
a different method, namely by evaluating the suitable 
Kubo formula for the spin current-spin current correla- 
tion function. By using the notation of Ref. d we have 
that 



(37) 



where j| and are the bare and dressed spin current 
vertices j| = {px/2m)T^ , J| = j| + F^. We then obtain 



cr 1 - (2apFr/uF)Tr(T=^r^) 



4e2 



1 + {2apFTY 



(38) 



The vertex corrections to the spin current vertex have 
been evaluated in Ref. (25j with the result F^ = 
vp{AapFT)~^T^ ■ One then obtains the vanishing of the 
spin current. 
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